Random field Ising model and community structure in complex networks 
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We propose a method to find out the community structure of a complex network. In this method 
the ground state problem of a ferromagnetic random field Ising model is considered on the network 
with the magnetic field B s = +oo, P>t = — oo, and Bi-t s ,t = for a node pair s and t. The ground 
state problem is equivalent to the so-called maximum flow problem, which can be solved exactly 
numerically with the help of a combinatorial optimization algorithm. The community structure 
is then identified from the ground state Ising spin domains for all pairs of s and t. Our method 
provides a criterion for the existence of the community structure, and is applicable to unweighted 
and weighted networks equally well. We demonstrate the performance of the method by applying it 
to the Barabasi- Albert network, Zachary karate club network, the scientific collaboration network, 
and the stock price correlation network. 
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I. INTRODUCTION 

The network theory is a useful tool for the study of 
complex systems. Universal features of some biologi- 
cal, social, and technological systems have been studied 
through their network structure 0, 0, El ■ Recent stud- 
ies have revealed that some complex networks have the 
community structure, which means that highly intercon- 
nected nodes are clustered in distinct parts. The com- 
munity may represent functional modules in biological 
networks 1, IE IE j industrial sectors in economic net- 
works [E a, and cliques of intimate individuals in social 
networks l(j ■ 

Recently various methods have been suggested for find- 
ing out the community structure in a given network . 
Girvan and Newman proposed an algorithm based on 
iterative removal of links with the highest betweenness 
centrality HE El 

The betweenness centrality of 
an edge is given by the number of the pathways pass- 
ing through it among shortest paths between all node 
pairs 13] . Nodes in different communities, if any, would 
be connected through rare inter-community links. Hence 
one could isolate communities by removing links with the 
highest betweenness centrality repeatedly. Similar meth- 
ods were also considered in Refs. [lJ,[lEllEl- Optimiza- 
tion techniques were also considered to find out the com- 
munity structure. In those approaches, the community 
structure is found by optimi zing an auxiliary quantity, 
such as the modularity 0, Il8j| . Some physical prob- 
lems turned out to be useful in detecting the community 
structure. For example, the g-state Potts model EH> the 
random walks j2jj, and the electric circuit problem pH 
were studied. 

Those methods proved to be successful in detecting ex- 
isting communities. On the other hand, it would be de- 
sirable to develop a method which can not only detect 
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the community structure but also verify its existence. 
Most algorithms developed are suitable for unweighted 
networks, whereas many real-world networks of interest 
are weighted . One may modify and generalize the al- 
gorithms developed for unweighted networks. However, 
such a generalization may not be straightforward |22| . So 
it would also be desirable to develop a method that works 
for unweighted and weighted networks equally well. 

In this paper we propose a method for finding out the 
community structure, which fulfills the requirements de- 
scribed above. Our approach is motivated from the ob- 
servation on the Zachary network, a classical example of 
social networks with the community structure |lfj. It 
is an acquaintance network of 34 members in a karate 
club. Once there arose a conflict between two influential 
members, which resulted in the breakup of the club into 
two. It is reasonable to think that the members would 
tend to minimize the number of broken ties, which can 
be accomplished by the breakup in accordance with the 
community structure. In fact, the resulting shape after 
the breakup coincides with the community structure of 
the original karate club network [h| . It suggests that the 
community structure of a given network may be found by 
simulating the breakup caused by an enforced frustration 
among nodes. 

We simulate the breakup by studying the ferromag- 
netic random field Ising model (FRFIM): The Ising spins 
(Tj = ±1 are assigned to all nodes i — 1,...,N, they in- 
teract ferromagnetically through links, and the quenched 
random magnetic field Bi is applied to each spin. The 
ferromagnetic interaction represents the cost for broken 
ties, and the random field is to introduce the frustration. 
In particular, we consider the case where the positive 
infinite magnetic field is applied to one spin and the neg- 
ative infinite magnetic field to another. It amounts to 
imposing the boundary condition that the two spins are 
in the opposite state. It simulates the conflict as raised 
by the two members in the Zachary network. Then, we 
will identify the community structure from the ground 
state spin domain pattern of the FRFIM. 
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This paper is organized as follows. In Sec. |H] we in- 
troduce the FRFIM in general weighted networks. The 
ground state problem of the FRFIM can be solved exactly 
with a numerical algorithm, which will be explained in 
Appendix. Then the method for finding out the com- 
munity structure is presented. In Sec. IIIII we apply 
the method to several networks and present the results. 
We conclude the paper with summary and discussion in 
Sec. El 



II. METHOD 

Consider a weighted network G of N nodes. Connec- 
tivity of G can be represented with the weight matrix 
{Jij\i,j = 1, • • • ,N}, where Jij is a prescribed weight 
or strength of a link between nodes i and j if they are 
connected or Jij — otherwise. We assume that the 
weights are non-negative, Jij > 0, and that the weights 
are symmetric, Jy = Jji. For an unweighted network, 
the matrix elements take the binary value or 1, and the 
weight matrix reduces to the usual adjacency matrix. 

The FRFIM on the network is defined with the Hamil- 
tonian 

h = ~^2 JijWj - Bi<Ti > w 

where a% = ±1 is the Ising spin variable at each node 
i. The spins interact ferromagnetically with the coupling 
strength {Jij}. They are also coupled with the quenched 
random magnetic field {Bi}. 

The FRFIM model has been studied extensively in d 
dimensional regular lattices in order to investigate the 
nature of the glass phase transition (see Ref. [23] and 
references therein). It has also been studied to investi- 
gate the disorder-driven roughening transition of inter- 
faces in disordered media [24| . The phase transition in 
the FRFIM on complex networks would also be interest- 
ing, which has not been studied so far. The issue will be 
studied elsewhere [25j . 

The specific feature of the FRFIM depends on the dis- 
tribution of the random field {Bi}. In this work, we 
consider the simple yet informative magnetic filed distri- 
bution given by 

{+00 , for i = s 
-00 , for i = t (2) 
, for i ^ s, t 

for certain two nodes s and t. It amounts to imposing the 
boundary condition that a s = +1 and a t = —1, which 
induces frustration among nodes. This specific random 
field distribution is adopted in order to mimic the con- 
flict as in the Zachary network. In the ground state nodes 
are separated into different spin domains, which will be 
related to the community structure of the underlying net- 
work. 




FIG. 1: Zachary karate club network. The links connecting 
nodes that are (not) in the same community are represented 
with solid (dashed) lines. The dotted lines separate the com- 
munities. 



As an explicit example, we consider the Zachary karate 
club network which is illustrated in Fig. ^ The node 
labeled as 1 (34) corresponds to the club instructor (ad- 
ministrator). They had a conflict, which resulted in the 
breakup. Nodes in the side of the administrator and the 
instructor after the breakup are denoted with circular 
and rectangular symbols, respectively. With Jy = 1 for 
all links and the magnetic field given by Eq. with 
s = 1 and t = 34, one can study the FRFIM on the net- 
work. Solving the ground state problem, we found that 
it has the degenerate ground states: The black (white) 
nodes belong to the + (— ) spin domain in all ground 
states, while the gray nodes (3, 10, 29) may belong to ei- 
ther domain. Note that the spin domains almost coincide 
with the actual shape of the network after the breakup; 
all black (white) nodes are in the side of the adminis- 
trator (instructor). The gray nodes are in the marginal 
state. It is reasonable to think that they do not belong 
to any community. In the previous work |Tol | , the node 3 
was misclassified. Our result hints that it is due to the 
marginality. 

The example clearly shows that the FRFIM is useful in 
finding out the community structure. For general appli- 
cation, (i) one needs to know the ground state(s) of the 
FRFIM of Eq. with the quenched random magnetic 
field given in Eq. for any node pair of s and t. Then, 
one needs to identify the set of all nodes that belong to 
the same spin domain as s and t in all ground states. 
Those sets will be called the cliques and denoted by C s 
and Ct, respectively. The number of nodes in the clique 
C will be called the clique size and denoted by \C\. (ii) 
More importantly, one needs to specify the node pair s 
and t which is relevant to the community structure. An 
arbitrary choice of s and t will not provide any informa- 
tion on the community structure. For example, if we take 
s = 12 and t = 15 in the Zachary network in Fig. ^ we 
obtain that C s = {12} and all other nodes are in Ct- This 
merely means that the node 12 is a peripheral node. 

For (i), the ground state problem of the FRFIM can be 
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solved exactly with the help of a numerical combinatorial 
optimization algorithm (see Appendix). This is achieved 
by mapping the ground state problem onto the minimum 
cut problem or the maximum flow problem p6|. The 
algorithm allows us to find all ground states, with which 
we can find the cliques C s and Ct for any pair of s and t. 
We explain the detailed procedure in Appendix. 

For (ii), the community structure can be found from 
the distribution of the clique sizes for all pairs of s and 
t. For a certain pair of s and t, one may have that 
\C S \ ~ \C t \ ~ 0(1) <C N. It happens when s and t are 
peripheral nodes of the network; most nodes are not in- 
fluenced by them. Such a pair does not provide any infor- 
mation on the community structure. One may have that 
\C S \ ~ 0(1) < \C t \ ~ 0{N). This happens when s is a 
peripheral node while t is inside the bulk. The cliques C s 
and Ct do not correspond to a community either. On the 
contrary, one may have that 0(1) <C \C S \ ~ \C t \ ~ 0{N). 
This happens only when there exist communities whose 
sizes are of the order of N, s and t are chosen among "in- 
fluential" nodes in different communities. In this case, 
we will regard the cliques C s and Ct as the communities 
in the network. 

In order to distinguish the different cases, we define the 
"separability" D s t for a node pair s and t as the product 
of the clique sizes, 



D. 



\C S \ ■ \C t 



(3) 



It ranges in the interval 1 < D st < N 2 /A. We propose 
that the community structure be detected with the dis- 
tribution of the separability D st for all pairs of s and t. 
If D st < 0{N) for all pairs of s and t, then we conclude 
that the network has no community structure. On the 
other hand, if D st ~ 0(N 2 ) for a certain pair of s and 
t, then we conclude that the network consists of commu- 
nities that can be identified from the cliques C s and Ct- 
Moreover, the nodes s and t may be regarded as influen- 
tial nodes of the communities. Therefore, in our method, 
the existence of the community structure is verified with 
the scaling behavior of the maximum value of the sepa- 
rability with the network size. 

For a given network size N, the scaling can be exam- 
ined with the quantity \nD st /\nN. Without the com- 
munity structure, it would be close to or much less than 
1 for all node pairs. A node pair with lnD st /\nN > 1 
indicates the presence of the community structure. 



III. RESULTS 



We test the method by applying it to the Barabasi- 
Albert (BA) network [22|, the Zachary karate club net- 
work [10|, the scientific collaboration network |lOj, and 
the stock price correlation network [2^ . In each network, 
the separability was calculated for all node pairs, and the 
separability distribution was examined with a so-called 
rank plot, where [In D/ In TV] is plotted against a normal- 
ized rank of each node pair. The rank is assigned to each 



( a ) 



c 



X 1 



( b ) 



0.5 1 

normalized rank 




0.5 1 

normalized rank 



( d ) 



J 



0.5 1 

normalized rank 



0.5 1 

normalized rank 



FIG. 2: The rank plot for the separability distribution for 
the BA network (a), the Zachary karate club network (b), 
the scientific collaboration network (c), and the stock price 
correlation network (d). 



node pair in the ascending order of the separability. It 
is then normalized so that the rank of the pair with the 
maximum value of the separability is equal to 1 . 

The BA network is an unweighted network. It is known 
that the BA network does not have a community struc- 
ture. We grew a BA network of N = 100 nodes, and 
calculated the separability D s t for all node pairs. The 
separability distribution is presented with the rank plot 
in Fig. |21 (a). We find that the separability is clus- 
tered at D st = 1 and near D st ~ N for all pairs of s 
and t, hence ]nD st / InN < 1. This confirms that the 
BA network does not have the community structure (see 
Fig. 121(a)), and demonstrates what the separability dis- 
tribution looks like for networks without the community 
structure. 

Next we study the separability distribution of the 
Zachary karate club network of N = 34 nodes, which 
is presented in Fig. (b). We found that D s t < N 
for all node pairs but (1,34) and (1,33). For the pairs 
(s, t) = (1, 33) and (1, 34), we obtained the same cliques, 
C s of 15 nodes and Ct of 16 nodes, which are marked 
with the black and the white symbols in Fig. ^ respec- 
tively (see also Fig. El (b)). Therefore we can conclude 
that there exist two communities in the network and that 
the node 1 is an influential node of one community and 
the nodes 33 and 34 are of the other community. In fact 
the nodes 1 and 34 correspond to the club instructor and 
the administrator, respectively. The detected communi- 
ties are in good agreement with the network shape after 
the breakup. 

We also investigate the community structure of a larger 
and more complex network. We examine the unweighted 
collaboration network of N = 118 scientists in the Santa 
Fe Institute 01 • l n tn i s network, two nodes (scientists) 
are linked if they coauthored at least one article. The 
rank plot is presented in Fig.|2|(c). One can see that the 
separability is distributed broadly, which indicates that 
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the network has multiple (more than two) communities. 

In such a case, the communities can be identified by 
applying our method hierarchically: First of all, one can 
find the node pair (so,io) with the largest separability, 
and the corresponding cliques C So and Ct . The clique 
may consist of a single community or be the union of 
several sub-communities. In order to investigate the sub- 
structure, one constructs the sub-network which consists 
of all nodes and links within each clique. Then, one can 
apply the method to the sub- networks. This can be per- 
formed hierarchically until a sub-network does not have 
the community structure any more. Or one may pro- 
ceed with the iteration only when the subnetwork size is 
equal to or larger than a threshold value m. The result- 
ing cliques can then be identified as communities up to a 
resolution m. 

With the hierarchical application of our method, we 
find the community structure of the scientific collabora- 
tion network as shown in Fig. [3] (c). Here, we identify 
all communities whose size are equal to or larger than 
m = 5. The community structure is in good agreement 
with that found in Ref . |l(J . 

Our method is also applicable to weighted networks. 
As an example of weighted networks, we study the eco- 
nomic network of 137 companies in the New York Stock 
Exchange market. The network is constructed through 
the stock price return correlation between the companies 
for the 21 year period from 1983 to 2003 28]. With the 
stock price Pi (t) of a company i at time t, the return is 
given by Ri(t) = lnP;(i + At) - \nPi(t) with the unit 
time interval At taken to be one day. Then, the stock 
price correlation is given by 

c _ ({Rj - {Ri))(Rj - (Rj))) 

J((Rj) - (Ri) 2 )((R?) {R 3 ?) ' 

where the angular bracket indicates the time average over 
the period. Its value ranges in the interval — 1 < CV, < 1, 
and is large for strongly correlated company pairs. It 
has been shown that the structural information of the 
economic system is encoded in the correlation matrix 

{(',,) BUS 

In order to apply our method, all weights are required 
to be non-negative. Hence, we assume that the weight is 
given by Jy = e aCij with a positive constant a taken to 
be 20. The weights are positive for all pairs of nodes, and 
the economic network is fully connected. The separabil- 
ity distribution is shown in the rank plot in Fig. El (d). 
As in the collaboration network, there are several non- 
trivial separability levels. We identified all communities 
whose size is equal to or larger than 3 with the same hi- 
erarchical method as in the collaboration network. The 
resulting shape of the network is illustrated in Fig. 01 (d). 
We confirmed that the communities are formed by com- 
panies in the same industrial sector. For example, the 
largest community consists of 13 companies in the energy 
sector. This study shows that our method works well 
for weighted networks. We note that many nodes (white 




FIG. 3: (Color online) The community structure of (a) the BA 
model network (N = 100), (b) the Zachary karate club net- 
work (N = 34), (c) the scientific collaboration network (N = 
118), and (d) the stock correlation network (N — 137). Nodes 
in different communities are distinguished with color. The 
white symbols represent the marginal nodes. 

symbols) remain unclassified. We attribute it to the fully- 
connectedness of the network. 



IV. CONCLUSIONS 

In this paper we have proposed the method for find- 
ing out the community structure of general networks. It 
is achieved by studying the ground state problem of the 
FRFIM on the networks with the magnetic field distri- 
bution given in Eq. J5J for two arbitrary nodes s and 
t. The cliques C s and Ct arc defined as the sets of all 
nodes that belong to the same spin domains as s and t in 
all possible degenerate ground states, respectively. The 
community structure is then manifested in the clique pat- 
tern for the pair with the maximum value of the separa- 
bility D st defined in Eq. ©. Our method is motivated 
from the observation on the Zachary karate club network, 
which shows that the resulting shape of the network af- 
ter breakup is determined by the underlying community 
structure. In our method, the response of the networks 
subject to schism is simulated with the FRFIM. 

In our method one can verify the existence of the com- 
munity structure of a given network with the scaling 
property of the separability: If the separability scales as 
D st < O(N) for all node pairs as in the BA network, the 
network does not have the community structure. On the 
other hand, if D st ~ 0(N 2 ) for a certain pair of s and 
t, one can conclude that the network has the community 
structure and that the nodes s and t are influential nodes 
in each community. Another advantage of our method is 
that it can be applied to both unweighted and weighted 
networks. Figure |3 shows the performance of the method 
in real- world networks. 

One of the weak points of our method is the time com- 
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plexity. Practically the ground state problem of the FR- 
FIM in sparse networks of N nodes has the time complex- 
ity of o\n 9 ) with ~ 1.2 .26]. Since one has to solve 
the ground state problems for all magnetic field distri- 
butions, the total time complexity scales as 0(N 2+ ). 
Hence, in the practical sense, our method is limited to 
networks of up to a few thousands of nodes. One may 
avoid the time complexity problem if the important nodes 
are known a priori. In the network theory, importance 
of nodes can be measured by, e.g., the degree or the be- 
tweenness centrality. Hopefully the community structure 
of large networks can be studied if one incorporates such 
importance measure into our method. 
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APPENDIX: MINIMUM CUT AND MAXIMUM 
FLOW PROBLEM 

This Appendix is intended to introduce the combina- 
torial optimization algorithm solving the ground state 
problem of the FRFIM. For more rigorous description, 
we refer the readers to Ref. ['2(i . 

Consider a network G of N nodes with the symmetric 
weight matrix {Jy > 0} — 1, • • • ,N). The ferro- 
magnetic random field Ising model on G is defined by 
the Hamiltonian in Eq. with the quenched random 
magnetic field {Bi}. The ground state is the spin con- 
figuration that has the minimum energy among all 2 N 
configurations. One might find the ground state by enu- 
merating all spin configurations, which is obviously time 
consuming and inefficient. We will explain the efficient 
way for solving the ground state problem. 

It is useful to introduce a capacitated network denoted 
by G': Having all nodes and links of G, G' contains two 
additional nodes S, called the source, and T, called the 
sink, and additional links between the source (sink) and 
the nodes with the positive (negative) magnetic field. G' 
is also a weighted network with the symmetric weight 
matrix {c Q( g} (a, (3 = S, T, 1, ■ • ■ , N). For a link (ij) from 
the original network G, the weight is given by cy = 2Jy. 
For the additional link the weight is given by csi = Bi for 
all i with Bi > and CjT = \Bi \ for all i with Bi < 0. The 
weight of the network G' is usually called the capacity. 
Figure 01 illustrates the relation between a network G of 
four nodes {a, b, c, d} and the corresponding capacitated 
network G'. 

In the capacitated network G' we define a ST-cut as a 
decomposition of all nodes into two disjoint sets S and T 
with S E S and T E T. It will be denoted by [S, T]. For 
a given [S, T], some links connect nodes in the different 
sets. The set of such links forms the boundary of the cut, 




(a) (b) (c) 

FIG. 4: (a) A network G with 4 nodes (filled circles) and 4 
links (lines) for the FRFIM. Figures represent the magnetic 
fields and the interaction strengths, respectively, (b) The cor- 
responding capacitated network G' with the link capacities, 
(c) The maximum-flow configuration with v* = 8. The dotted 
lines represent saturated links with x* a g — c a p. The dashed 
lines X, Y , and Z represent boundaries associated with ST- 
cut s. 

which is denoted by (S, T) = {(a(3)\a E S,j3 E T}. The 
cut capacity C[S, T] is then defined as the total sum of 
the capacity of the boundary links, that is, 

C[S,T\= Yl c ^ ' ( A - : ) 

(a/3)e(S,T) 

Figure 0] shows some examples of the cut. The boundary 
denoted by X is associated with a cut [{S, a, b, c}, {T, d}] , 
whose cut capacity is 14. 

There exists one-to-one correspondence between the 
Ising spin configuration on the weighted network G and 
the cut [S,T] of the capacitated network G'. It is 
achieved by assigning Ui = +1(— 1) for all nodes i in S(T) 
and vice versa. Hence, the sets S and T correspond to 
up and down spin domains, respectively, and the bound- 
ary (<S, T) corresponds to the spin domain wall. Fur- 
thermore, one can easily verify that the energy E of the 
FRFIM of a spin configuration {a;} and the cut capacity 
C[S, T] satisfy the relation 

E({a i }) = C[S ) T]+E (A.2) 

where E = - £\ ^ Jy/2 - J2i I- 8 *!/ 2 - Therefore, solving 
the ground state of the FRFIM on G is equivalent to 
finding the optimal ST-cut on G' whose cut capacity is 
minimum. It is called the minimum cut problem. 

The minimum cut problem can be further mapped on 
to the maximum flow problem: On the capacitated net- 
work G", a flow is to denote a set of flow variables {cc a( a} 
defined for all links in G' which are subject to a capacity 
constraint 

< x af3 < c a p (A. 3) 

and a mass balance constraint 

2^ x a /3 - xp a = v5(a, S) - v6(a,T). (A.4) 

Here J^' means a sum over all adjacent nodes of a, S() 
denotes the Kronecker 5 symbol, and v is a non-negative 
parameter. The mass balance constraint allows us to 
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interpret the flow {i^} as a conserved flux configuration 
of a, e.g., fluid which is originated from the source S by 
the amount of v and targeted to the sink T through the 
network G' . 

Due to the capacity constraint, there exists the upper 
bound in v, beyond which a flow satisfying Eqs. I|A.3(I 
and (|A.4|I does not exist. Then, the question that arises 
naturally is to find the maximum value v* and the cor- 
responding flow {x* a p\ that can be delivered. This is the 
maximum flow problem. 

The celebrated max-flow/min-cut theorem of Ford and 
Fulkerson [3(j states that for a given capacitated network 
G", the maximum flow v* is equal to the minimum cut 
capacity, that is to say, 

v* = min CIS, T] . (A.5) 

[S.T] 

The rigorous proof of the theorem can be found else- 
where 26] . Intuitively the theorem states that the max- 
imum flow is limited by the bottleneck in the network 
whose capacity is given by the minimum cut capacity. 

The maximum flow problem can be solved numerically 
in a polynomial time with the augmenting path algo- 
rithm or the preflow-push/relabel algorithm j2|||3(]]. In 
the augmenting path algorithm, one repeatedly searches 
for a path from S to T via unsaturated (x a p < c a p) links 
and updates {x a p} by augmenting flows along the path. 
When the augmenting path does not exist any more, the 
resulting flow corresponds to the maximum flow config- 
uration. The preflow-push/relabel algorithm is a more 
sophisticated and efficient algorithm. 

Once the maximum flow configuration {x* a a} is found, 
the minimum cut is constructed easily. Let Ss be the set 
of all nodes of G' that can be reachable from the source 
S only through unsaturated (x* a p < c a p) links. Trivially, 
Ss does not include the sink T, since there does not exist 
any augmenting path in the maximum flow configuration. 
Hence, the set Ss and its complement Ss defines a cut 
[5s, 5s], which is indeed a minimum cut of G' . 

One may find the minimum cut alternatively. Let Tt 
be the set of all nodes of G' that can be reachable from 
the sink T only through unsaturated links. Then, Tt and 
its complement Tt defines a cut [Tp^Tr], which is also 
the minimum cut. 



The two cuts [Ss, Ss] and [Tt,Tt] may be different, 
which implies that the corresponding FRFIM has degen- 
erate ground states. In that case, all degenerate ground 
states can be found systematically |26| . In this work, we 
are interested in the spins that are fixed in all ground 
states. One can easily verify that all nodes i £ Ss (Tt) 
except for S (T) are in the spin state cr, = +1 (— 1) in all 
ground states. The other nodes j Ss and Tt may be 
in either state <jj = ±1. 

We provide an example illustrating the mapping be- 
tween the FRFIM and the maximum flow or the min- 
imum cut problem in Fig. 0] The maximum flow con- 
figuration is depicted in Fig. 0] (c) with the maximum 
flow v* = 8. The links drawn with dotted lines are satu- 
rated (x* a p = c a p). The sets of all nodes that are reach- 
able from S and T through unsaturated links are given 
by Ss = {S, a} and Tp = {T, &, d}. They yield the mini- 
mum cuts [iSs,iSs] and [Tt,Tt] whose boundaries are Y 
and Z, respectively. Hence, one finds that a a = +1 and 
<?b = (Td = — 1 in all degenerate ground states. The node 
c does not belong to neither Ss nor Tt- Hence er c may 
be either +1 or —1. 

In the present work, we consider the FRFIM on a 
weighted network G with the specific magnetic field dis- 
tribution given in Eq. J2Jl for a certain node pair s and 
t. Then, we need to find the clique C s (fit) of s (t) which 
is the set of all nodes that are in the same spin state as 
s (t) in the ground state. We summarize the method to 
find the cliques: 

1. Construct the capacitated network G' . 

2. Find the maximum flow configuration {x*^} using 
the numerical algorithms. 

3. Find the set Ss (Tt) of all nodes that are reachable 
from S (T) through unsaturated links with x* a ^ < 

Cap- 

4. Then, the cliques are given by C s = Ss — {S} and 
C t = T T -{T}. 

After finding the cliques, the community structure can 
be investigated with the method explained in Sec. [n] 
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